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Abstract 

We carry out a high-precision Monte Carlo simulation of the two-dimensional 
0(3)-invariant ci-model at correlation lengths .^ up to ~ 10^. Our work employs 
a new and powerful method for extrapolating finite- volume Monte Carlo data 
to infinite volume, based on finite-size-scaling theory. We discuss carefully the 
systematic and statistical errors in this extrapolation. We then compare the ex- 
trapolated data to the renormalization-group predictions. The deviation from 
asymptotic scaling, which is ~ 25% at .^ ~ 10^, decreases to ~ 4% at .^ ~ 10^. 

PACS number(s): ll.lO.Gh, 11. 15. Ha, 12.38.Gc, 05.70.Jk 



Two-dimensional nonlinear cr-models are important "toy models" in elementary- 
particle physics because they share with four- dimensional nonabelian gauge theories 
the property of perturbative asymptotic freedom [1]. However, the nonperturba- 
tive validity of asymptotic freedom has been questioned [2]; and numerical tests of 
asymptotic scaling in the 0(3) cr-model at correlation lengths i^ ~ 100 have shown 
discrepancies of order 25% [3,4]. In this Letter we employ a new finite-size-scaling 
extrapolation method [5] (see also Liischer et al. [6] and Kim [7] for related work 
[8]) to obtain high-precision estimates (errors < 2%) in the 0(3) cr-model at corre- 
lation lengths 1^ up to ~10^. We find that the discrepancy has decreased to sa 4%, 
in good agreement with the asymptotic-freedom predictions. 

We study the lattice cr-model taking values in the unit sphere S^~^ C R , with 
nearest-neighbor action 'H{(t) = —jSY^f^x ■ f^y Perturbative renormalization-group 
computations predict that the (infinite-volume) correlation lengths ^y-^P) and i^^^^ 
[9] beh. 
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as /3 — )• oo. Three-loop perturbation theory yields [11,12] 

a, = -0.014127+ Q-^)/(iV- 2). 
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The nonperturbative constant (7^(exp) has been computed recently using the ther- 
modynamic Bethe Ansatz [13]: 
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The remaining nonperturbative constant is known analytically only at large N [14]: 
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Previous Monte Carlo studies up to i^ ~ 100 agree with these predictions to within 
about 20-25% for iV = 3, 6% for iV = 4 and 2% for iV = 8 [4,15]. 

Our extrapolation method [5] is based on the finite-size-scaling Ansatz 
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where O is any long-distance observable, s is a fixed scale factor (usually s = 2), L 
is the linear lattice size, Fq is a universal function, and cu is a correction-to-scaling 
exponent. We make Monte Carlo runs at numerous pairs {fi^L) and (jS^sL); we 
then plot C(/3, sL)/0(l3^ L) versus i^(/3, L)/ L^ using those points satisfying both 
i^(/3,-L) > some value ^min and L > some value Lmin- If all these points fall with 
good accuracy on a single curve, we choose a smooth fitting function Fq- Then, 
using the functions F(_ and Fq^ we extrapolate the pair (C^O) successively from 



L ^ sL ^ s^L —)•...—)• oo. See [5] for how to calculate statistical error bars on 
the extrapolated values. 

We have chosen to use functions Fo of the form 

Fo{x) = 1 + aie-i/" + a^e-^l'' + . . . + a^e"-^/" . (6) 

This form is partially motivated by theory, which tells us that F{x) — )• 1 expo- 
nentially fast as x — )• [16]. Typically a fit of order 3 < n < 12 is sufficient; we 
increase n until the x^ of the ffi becomes essentially constant. The resulting x^ 
value provides a check on the systematic errors arising from corrections to scaling 
and/or from inadequacies of the form (6). The discrepancies between the extrap- 
olated values from different L at the same /3 can also be subjected to a x^ test. 
Further details on the method can be found in [5]. 

We simulated the two-dimensional 0(3) cr-model, using the Wolff embedding al- 
gorithm with standard Swendsen- Wang updates [17,18,10]; critical slowing- down ap- 
pears to be completely eliminated. We ran on lattices L = 32, 48, 64, 96, 128, 192, 256, 
384,512 at 180 different pairs {fi^L) in the range 1.65 < /3 < 3.00 (corresponding 
to 20 ^ i^oo ^ 10^). Each run was between 10^ and 5 X 10^ iterations, and the total 
CPU time was 7 years on an IBM RS-6000/370. The raw data will appear in [19]. 

Our data cover the range 0.15 < ^(L)/L < 1.0, and we found tentatively that 
a tenth-order ffi (6) is indicated: see Table 1. Next we took ^min = 20 and 
sought to choose Lmin to avoid any detectable systematic error from corrections 
to scaling. There appear to be weak corrections to scaling (< 1.5%) in the region 
0.3 ^ ^(L)/L ^ 0.7 for lattices with L ^ 64-96: see the deviations plotted in 
Figure 1. We therefore investigated systematically the x^ of the ffis, allowing a 
different Lmin for ^(L)/L < 0.7 and > 0.7: see Table 1. A reasonable x^ is obtained 
when n > 9 and Lmin ^ (128, 64). Our preferred ffi is n = 10 and Lmin = (128, 64): 
see Figure 2, where we compare also with the perturbative prediction 
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vahd for x > 1, where a = l/(iV - 1), wq = {N - 2)/27r and wi = {N - 2)/(27r)2. 

The extrapolated values i^^^^ from different lattice sizes at the same /3 are con- 
sistent within statistical errors: only one of the 24 fi values has a x^ too large at 
the 5% level; and summing all fi values we have x^ = 86.56 (106 DF, level = 92%). 

In Table 2 we show the extrapolated values i^^^^ from our preferred fit and some 
alternative ffis. The discrepancies between these values (if larger than the statistical 
errors) can serve as a rough estimate of the remaining systematic errors due to 
corrections to scaling. The statistical errors in our preferred ffi are of order 0.2% 
(resp. 0.7%, 1.1%, 1.6%) at 1^00 ~ 10^ (resp. 10^, 10"*, 10^), and the systematic errors 
are of the same order or smaller. The statistical errors at different fi are strongly 
positively correlated. 

In Figure 3 (points -|- and X ) we plot Coo\,gtiynate (128 64) divided by the two-loop and 
three-loop predictions (l)-(4). The discrepancy from three-loop asymptotic scaling, 
which is ^ 16% at /3 = 2.0 (^ ~ 200), decreases to ^ 4% at /3 = 3.0 (^ ~ 10^). 



This is roughly consistent with the expected l//3^ corrections. The shght bump at 
2.3 ^ /3 ^ 2.6 is probably spurious, arising from correlated statistical or systematic 
errors. 

We can also try an "improved expansion parameter" [20,4,12,19] based on the 
energy E = {(Tq ■ <Ti). First we invert the perturbative expansion [21,12] 






N-1 N-1 0.005993(iV- 1)2 + 0. 007270(iV- 1) 
1/3 32/32 Js 

(8) 
and substitute into (1); this gives a prediction for i^ as a function of 1 — £^. For E 
we use the value measured on the largest lattice; the statistical errors and finite-size 
corrections on E are less than 5 X 10~^, and therefore induce a negligible error (less 
than 0.5%) on the predicted i^. The corresponding observed/predicted ratios are 
also shown in Figure 3 (points □ and O). The "improved" 3-loop prediction is in 
excellent agreement with the data. 

Let us summarize the conceptual basis of our analysis. The main assumption 
is that if the Ansatz (5) with a given function F(_ is well satisfied by our data for 
Lfnin 1^ L < 256 and 1.65 < /3 < 3, then it will continue to be well satisfied for 
L > 256 and for /3 > 3. Obviously this assumption could fail, e.g. if [2] at some 
large correlation length (> 10"^) the model crosses over to a new universality class 
associated with a finite-/? critical point. In this respect our work is subject to the 
same caveats as any other Monte Carlo work on a finite lattice. However, it should 
be emphasized that our approach does not assume asymptotic scaling [eq. (1)], as 
/3 plays no role in our extrapolation method. Thus, we can make an unbiased test 
of asymptotic scaling. The fact that we confirm (1) with the correct nonperturbative 
constant (3)/(4) is, we believe, good evidence in favor of the asymptotic-freedom 
picture. We are unable to imagine how, if there were in fact a finite-/? critical 
point [2], the "preasymptotic" region at /3 < 3 would mimic not only asymptotic 
freedom but also the nonperturbative constant predicted by the thermodynamic 
Bethe Ansatz. 

Details of this work, including an analysis of the susceptibility x, will appear 
elsewhere [19]. 
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ing us to push to ever larger values of ^/ L. We also thank Steffen Meyer, Adrian 
Patrascioiu, Erhard Seller and UUi Wolff for helpful discussions. The authors' re- 
search was supported by CNR, INFN, DOE contracts DE-FG05-85ER250000 and 
DE-FG05-92ER40742, NSF grant DMS-9200719, and NATO CRG 910251. 
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Figure 1: Deviation of points from fit to F(_ with s = 2, ^min = 20, Lmin = 128, 
n = 10. Symbols indicate i: = 32 (+), 48 (*), 64 (x), 96 (?8^), 128 (□), 192 (^), 256 
(O). Error bars are one standard deviation. Curves near zero indicate statistical 
error bars (it one standard deviation) on the function F^(x). 
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Figure 2: ^(/3, 2i:)/^(/3, L) versus C{J3,L)/L. Symbols indicate L = 32 (+), 48 (*), 
64 (x), 96 {^), 128 (□), 192 ()p(), 256 (O). Error bars are one standard deviation. 
Solid curve is a tenth-order fit in (6), with ^min = 20 and Lmin = 128 (resp. 64) for 
^{L)l L < 0.7 (resp. > 0.7). Dashed curve is the perturbative prediction (7). 
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Figure 3: £ ,• , noc c.l^/£L +(,.™ versus /3. Error bars are one standard deviation 

o ^oo, estimate (128,64)/ ^oOjt/ieor / 

(statistical error only). There are four versions of C^o t/ieor- standard perturbation 
theory in 1//3 gives points + (2-loop) and X (3-loop); "improved" perturbation 
theory in 1 — £^ gives points □ (2-loop) and O (3-loop). 
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Table 1: x^ and confidence level for the fit (6) of i^(/3, 2L)/i^(/3, L) versus i^(/3, L)/ L. 
DF = number of degrees of freedom. The first (resp. second) Lmin value applies for 
C{L)/L < 0.7 (resp. > 0.7). In aU cases ^rmn = 20. 
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1402 ( 22) 


2499 ( 41) 
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— 


— 


— 





Table 2: Estimated correlation lengths i^^^^ as a function of /3, from various ex- 
trapolations. Error bar is one standard deviation (statistical errors only). All 
extrapolations use s = 2, ^min = 20 and n = 10. The first (resp. second) Lmin value 
applies for ^(L)/L < 0.7 (resp. > 0.7). Our preferred fit is Lmin = (128,64), shown 
in italics. Kim is the estimate from [7]. 
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